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Abstract We use Permi coordinates to calculate the canonical partition func- 
tion for an ideal gas in a circular geodesic orbit in Schwarzschild spacetime. 
To test the validity of the results we prove theorems for limiting cases. We 
recover the Newtonian gas law subject only to tidal forces in the Newtonian 
limit. Additionally we recover the special relativistic gas law as the radius of 
the orbit increases to infinity. We also discuss how the method can be extended 
to the non ideal gas case. 
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1 Introduction 

Contributions to the development of a theory of relativistic Statistical mechanics, in- 
cluding attempts to find a fully covariant theory, date back to the early part of the 20th 
Century. A useful summary of contributions up to 1967 is given by Hakim Hl. Since 
that time distinct proposals for a general framework for relativistic Statistical mechan- 
ics have continued to emerge. Por example, Horwitz, Schieve and Piron |2 1 constructed 
classical and quantum Gibbs ensembles in the context of special relativity using "his- 
torical time" and a statistics of events off the mass shell. Miller and Karsch |3| applied 
the constraint formalism to many-particle special relativistic Systems in Statistical me- 
chanics in Order to include two-particle interactions. Droz-Vincent [4J considered a 
statistics of orbits consistent with Hakim's approach including a tentative definition of 
equilibrium. More recently, Montesinos and Rovelli Q proposed a theory of Statistical 
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mechanics that is not necessarily tied to the notion of energy or a preferred time axis 
and which is consequently not necessarily governed by the notion of temperature. 

In this paper, we consider an ideal gas enclosed in a Container in orbit around the 
central mass in Schwarzschild spacetime. It is assumed that the gas particles, along 
with the Container, do not significantly alter the background Schwarzschild metric.We 
develop a Statistical mechanical formalism that yields approximations of the canoni- 
cal Helmholtz free energy for the volume of ideal gas. These approximations involve 
cross terms between space coordinates, momentum coordinates, and Riemann curva- 
ture tensor components, but the expressions are still computationally accessible enough 
to allow comparision to their special relativistic and Newtonian analogs. We test the 
validity of these expressions by proving rigorous results for special relativistic and 
Newtonian limits. Technical issues related to the existence of the thermodynamic limit 
are discussed in Remark 1 below. 

Our method begins with the calculation of the Schwarzschild metric in Permi coor- 
dinates along the circular geodesic orbit. For that purpose, we use the results of |7 1 and 
||8|. The choice of a circular orbit in the spherically Symmetrie Schwarzschild space- 
time makes it plausible on physical grounds that the gas will achieve equilibrium. Since 
Permi coordinates are Minkowskian to first order near the orbit, the Statistical mechan- 
ical partition function should closely resemble the analogous well-known expression 
in special relativity (see e.g. |9|), but with corrections to second order We note that 
the term "ideal gas" is somewhat misleading in the context of general relativity. This is 
because a volume of gas subject to no forces is still affected by the curvature of space- 
time, and this corresponds to a Newtonian gas subject to tidal forces, but otherwise 
"ideal." We elaborate on this correspondence in See. 4. 

Our formulation and results may be compared with |5 1 and |6|. In Section 5 of O 
it was shown that even for the ideal gas in Minkowski spacetime, "there is no Lorentz- 
invariant thermal State: a thermal State is in equilibrium in a preferred Lorentz frame 
of reference, and therefore breaks Lorentz invariance." It was further argued in [5 1 that 
when the gas is enclosed in a box, the coordinate System in which the center of mass 
of the System is at rest is the preferred Lorentz frame, and the canonical distribution 
density for an ideal gas in Minkowski spacetime then reduces to a product of terms 
of the form expßpo (where po is the negative of the energy of a particle and ß is the 
inverse temperature). Analogously, we work in the locally Minkowskian Permi normal 
coordinate System of the rest frame of the box of gas in Schwarzschild spacetime. We 
expect this to be the preferred reference frame for equilibrium. In Section 8 of |6], 
Chernikov gave a Solution to the kinetic Boltzmann equation for the single-particle 
equilibrium distribution of an ideal gas in the Schwarzschild coordinates. In contrast 
to the methods of [6J, we use Statistical mechanical arguments for derivation of the 
canonical partition function. We make further comparisons in Section 5. The time- 
dependence of the Permi metric considered here necessarily introduces complications. 

In See. 2 we introduce notation and compute the Schwarzschild metric in Permi 
coordinates to second order. In See. 3 we derive the volume form for phase space, give 
a theorem for coordinate independence, and establish a local Version of Liouville's 
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Theorem on a 6A^-dimensional phase space for particles. We then find approxima- 
tions for the partition function and free energy for the ideal gas. The proofs of the 
main resuhs of this section appear in Appendices A and B. In See. 4 we prove that 
the non relativistic Newtonian limit of our general relativistic free energy function is 
the expected result corresponding to classical tidal forces. We also prove that the limit 
of our general relativistic estimate of the free energy, as the radius of the orbit goes to 
infinity, equals the free energy for the special relativistic case found by Pauli 121 . See. 
5 is devoted to conclusions. 



2 Fermi coordinates 

Throughout we use the sign Conventions of Misner, Thorne and Wheeler IfTÜI and ex- 
press the Schwarzschild line element in terms of the usual variables as, 

ds^ = -(l- dt^ + ^ + r^ide^ + sin^ Od^^). (1) 

r 

Consider a Container of gas molecules in orbit around the central mass M in Schwar- 
zschild spacetime. We assume that neither the gas nor the Container significantly con- 
tribute to the gravitational field, so that the Schwarzschild metric is effectively unal- 
tered by their presence. We assume further that a center of mass of the Container of the 
gas may be identified, and that it follows (or very nearly follows) a circular geodesic 
orbit ^ of radius r > with = 7i;/2; (a discussion of the existence of a center of 
mass related to our circumstances is given in ifTTl . see in particular p. 209). The four- 
velocity tangent vector along is {i,Q,Q,id(l) /dt). In a neighborhood of the circular 
orbit all spacetime points may be charted using Fermi coordinates whose origin 
is the center of mass of the Container of the gas. Parker and Pimentel [[8| calculated 
the Fermi normal basis vectors, or tetrad, for a circular geodesic orbit at radius r and 
gave some of the Riemann tensor components in that tetrad. Below, we give the Fermi 
basis vectors (after correcting a misprint in Eq. (5.24) of lH)) and all of the nonvanish- 
ing Riemann tensor components in the Fermi frame. Hereafter barred indices refer to 
the Schwarzschild coordinates of ([T]) while unbarred indices refer to Fermi coordinates 
{x^ ,x^}. As usual Greek indices run over 0, 1,2,3 and Latin over 1,2,3. Let, 

/ ( f fl i4n / <5r dr dO d(j) 
= = (2) 

Then a parallel tetrad for Fermi coordinates is. 



^Since circular orbits are stable only for r > 6M, it may be argued that equilibrium is possible only for 
such orbits. 



3 



e / 



eo = (-,0,0,-), (3) 



^1 = 1-^:^'^ ^ (4) 

62 = (o,0,i,0^, (5) 
' " i rXV2 sin(a(/)),0, ^^^^^ I , (6) 



where. 



2M /r-3M 

^=1 , a = J , (7) 

r \ r 

e is the energy per unit mass of a test particle following the orbit of the center of mass 
of the Container of gas, and l is its angular momentum per unit mass. Thus, 



r-2M / M 

£ = , , and l = r\ — . (8) 

y/r{r-3M) \ r-3M 

The tetrad vectors may be expressed expHcitly in terms of the proper time = T of the 
Permi observer. On the circular geodesic "if we have that 

<^ = 4' (9) 

where the dot indicates the derivative with respect to the proper time along this geodesic, 
namely T. For the geodesic , l is given by ([8]), therefore it foUows that 



^ = (10) 

We now use the relation 

i^aßrS=RKlßv'^Seßel'e^, (11) 
to obtain the Riemann tensor components in the Permi basis (cf. fSl). 
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M{r + 3{r-2M)cos{2ac^)) 
2(r-3M)r3 ' ^^^^ 



3M{r-2M) sin(2a0) 
2{r-3M)r^ 



^0103 = ^TT TTTT^ ' (1^) 



M(r-3(r-2M)cos(2a^)) 
2(r-3M)r3 



^0303 = TT ^TTT^ 1 (15) 



3M3/2(r-2M)>/2cos(a^) 

^«^'^ - (r-3M)r3 ' ^^^^ 

3M3/2(r-2M)i/2 sin(a^) 
"""^^^ ^ ^-3Mp ■ ^^^^ 

From (|7]l and (fTÜb we have that 

a0 = ^T. (18) 

(fTSl l expresses the argument of the trigonometric functions appearing in the the curva- 
ture components, (fT2] |-(fT7l). in terms of T. 

In addition to the usual symmetries which connect the various components of the 
Riemann tensor we find the following relations 

^0313 = -^0212, ^0223 = — ^0113j (19) 
^1212 = -^0303. ^1313 = -^0202, (20) 
^1223 = -^0103: ^2323 = -^0101- (21) 

The metric components in Permi coordinates to second order are given (l?], ifTÜl . 
and mi) by. 
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^00 = -\-RomJr + 0{x^), (22) 

^0,- = /' = -|7;o,™xV" + 0(x3), (23) 

gij - 5,j-^-RujJ^ + 0{x^), (24) 

= -l+Äo/OmxV« + 0(x3), (25) 
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= 5" + -Rnjmx'x"' + 0{x'), (26) 
^ = -l + ^{Rim~2Rmm)x'x'" + 0{x^), (27) 

whereg = det(§„ß). 

3 Phase space and the canonical partition function 

The Hamiltonian for a System of ideal gas particles is given by 

N In 

where H, ^ -gf p,aP,ß, (28) 



Here H, is the Hamiltonian for the /th particle. For notational simplicity we shall omit 
the subscript / labeling the particle coordinates and momenta, when it is clear that we 
are referring to a one-particle Hamiltonian. The State of a Single particle consists of its 
four spacetime coordinates together with its four-momentum coordinates {x"' ^pß}. If 
the proper mass, m > 0, of the particle is fixed, then the one-particle Hamiltonian is 
given by, 

H, = \8"PpaPß^-Y- (29) 

Eq. ( |29] l allows us to reduce the dimension of the State space to seven. The energy 
of a test particle with four-momentum pß according to an observer with four-velocity 
u" in Schwarzschild spacetime, is —u"pa, which is a scalar invariant under coordinate 
transformations. In the Permi coordinates described in the previous section, an observer 
with fixed space coordinates {x^} and time coordinate T, e.g., the observer at any fixed 
location in the Container of the orbiting gas, with four-velocity m" — (1,0,0,0), will 
measure energy — po of a particle with four-momentum pß at that point. With the 
assumption po < we obtain. 
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-g^'p. + Vis"' 8"' - 8V)PjPk - ^ 
/?o = • (30) 

Let P be the sub bündle of the cotangent bündle of spacetime with each three di- 
mensional Aber determined by (l3Ül i. 

Proposition 1 The volume 7-form cb on P given by, 

Ö) = \dx^ Adx^ Adx^ Adx^ Adpi Adp2Adp3, (31) 
P 

is invariant under all coordinate transformations. 

The proof of Proposition I which is analogous to a result in lfT3l (see also lfT4j . p. 264 
for generalizations) is given in Appendix A. 

We next restrict our attention to the Permi Observer along the circular geodesic 
with fixed T. Phase space for a Single particle is determined by the space slice orthogo- 
nal to the Permi Observer's four-velocity, with the associated momentum coordinates, 
along with the volume form given by the interior product (cf. fT?!, p. 3) i{d /dT)ä) of 
the vector (9/(9t with ö). Then, 

mi{d/dx)Ö} = dx^ Adx^ Adx^ Adpi Adp2Adp3 

12 3 

= dx dx dx dpidp2dp3. (32) 

It follows from Proposition 1 that the 6-form given by ( |32] | is invariant under coordinate 
changes of the space variables (only), with t fixed. Physically this means that the 
calculations that follow below are independent of the choice of the Permi tetrad, (|2]i- 
(|6]l. We note in particular that with T fixed, po is invariant under changes of the space 
coordinates. 

Liouville's theorem for the invariance of phase space volume under the dynamics of 
particle motion is one of the foundations of Statistical mechanics. In order to establish 
a local Version of Liouville's theorem associated with a product of 6-forms each of the 
form given by ( [32] i we return briefly to the cotangent bündle of spacetime and consider 
the symplectic form, röo, given by, 

cbo = dx^^ Adpo + dx^ A dpi +dx^ Adp2 + dx^ Adpj. (33) 

The Hamiltonian, Hj, governing the dynamics is given by (|29] l. Henceforth, for the 
duration ofthis article, we use gfj,ß to represent the second order approximation to the 
Permi metric given by \22^ - \24^ rather than the exact ( uncalculated) values. Likewise, 
we use only the second order approximations for the components of the inverse metric 
andg"'^ . 
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The metric and hence H depend on only through the Riemann tensor compo- 
nents. It follows from (fT2]l-(fT7]i that the Hamiltonian is given to second order in all 
four spacetime coordinates {x^ ,x^,x^} by setting x" = 0. In fact, expanding to third 
Order in the spacetime coordinates, it is not difficult to show that the time derivatives 
of the metric tensor Clements are much smaller than the derivatives in the spatial di- 
rections at spacetime points close to the geodesic circular orbit (cf. |7|, Section IX). 
For the purpose of approximating Statistical mechanical quantities we therefore make 
a further approximation and set x*^ = in (fT2b-(fT7]i. This amounts to selecting a point 
on the circular geodesic orbit and assuming an approximately time independent Hamil- 
tonian. In light of the comments following Eq. ( |32l ). we would obtain the same final 
results from the selection of any other point on the orbit, and we briefly discuss this in 
Section 4 below. We therefore identify po = —E < 0, where the energy £ of a particle 
is positive and constant, and it follows that we can reduce the dimension of phase space 
by two by eliminating po and its conjugate cyclic coordinate.jc". 

We next consider an ideal gas consisting of particles in a one-dimensional box; 
the generalization to additional degrees of freedom is straightforward and we omit it in 
Order to keep our indices simple. We label the coordinate, momentum, and proper time 
of the Ith particle, (/ = 1, . . . ,A^) by qi,pi, T/, respectively. We assume that T, ~ T;(t), 
where T the proper time of the center of mass of the System (in the interval between 
coUisions). Note that the function T/(t) may be any smooth function of T. 

We may dehne the Hamiltonian vector field, vh, as follows 



VH 



( ^1 \ 
P\ 

V PN j 



dx 

dpi 
dl 



dqN 
dT 

dpN 
dT 



/ dqi dX\ \ 
dxx dx 

dpi dTi 
dTi dx 



dqN dTM 
dXN dx 

dpN dXN 
\ dXN dx / 



(34) 



Liouville's theorem, which states that the Hamiltonian System is measure-preserving, 
follows from ( l35] l: 



div vh = 



^^^P]_^ I f^qN 

dqi dpi dqN 
d f dqA 
dqi \dXi 

dqidpi dpidqi 



dp 



N 



dp 



N 



d 
dpi 
d^H 



dpi 
dxi 



dXi 
dx 



dXN 

'dx' 



dxi 




dXN 


— + ••• + 
dx 




dx 



= 0. 



(35) 
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We thus have a local Version Liouville's theorem for our gas restricted to a small vol- 
ume in circular orbit, that is, we may assert that Liouville's theorem holds to low order 
in the spacetime variables. 

In Order to proceed to the next step and dehne the thermodynamic free energy cor- 
responding to g^ß (see (l38l l below) we will need the following theorem. The proof is 
given in Appendix B. 

Theorem 1 There exists a positive function f{M,r) definedfor r > 3M > such that: 

(a) f{M,r) isadecreasingfunctionofM and an increasing function ofr and f{M,r) — > 
ooasr—foaorM—i-0 when the other variable is heldfixed. 

(b) If\x'\ <f{M,r) fori ^1,2,3, then 

Po<-fe(bil + b2l + b3l), (36) 

for sotne positive number k, where po is the function of space and momentum coordi- 
nates defined by diOt . 



The canonical partition function for the gas under consideration depends on temper- 
ature. One would expect small variations in the temperature T within the gas volume 
(a) because of pressure variations due to gravity and (b) because of tidal effects on 
the particles in a Container in orbit around M (cf. ifTSl . ifTöl ). These effects, although 
very small, would also be present in the nonrelativistic case. We assume that our Sys- 
tem is in contact with a 'heat bath' and consequently in thermal equilibrium to a good 
approximation. Following the usual Convention, we write, 

ß-^, (37) 

where T is the temperature of the gas in volume V and k is Boltzmann's constant. 
Turning now to the canonical partition function, we write, 

dx = dx^dx^dx^ , dp — dp\dp2dp^, 
dx'^dp'^ = dx\dx\dx\dp\idpi2dp\-i ■ ■ ■dxlfdxljdxl/dpMidpN2dpM3- (39) 

where is the number of particles, and {po)i is given by (|3Ü] l with the coordinates and 
momenta corresponding to the Ith particle. The volume integrals in (1381 ) have limits 
of integration determined by the ränge of Permi coordinates that defines the volume of 
the gas. 

The following result follows immediately from Theorem 1 . 
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Corollary 1 The partition function given by \38i is finite for any particle number N, 
any inverse tempemture ß, and any volume V for which \x' \ < f{M,r) for i— 1,2,3, 
where f{M, r) is given by Theorem 1. 

Remark 1 The physical significance of the bound on |;ic'| given by f{M,r), a function 
that increases with increasing r or decreasing M, is that the Permi coordinates of the 
volume of the Container of the gas may be increased - and therefore the volume itself 
may be enlarged - as the orbital radius increases, or the mass of the central body de- 
creases. Informally, in the limit as r ^ oo orM — > 0, the volume of the gas may increase 
to infinity, and thus the thermodynamic limit may be considered with the bounds ( [36] l 
on the integrand of the partition function ( |38] |. We leave for future investigations the 
question of whether or not Theorem 1 can be sufficiently sharpened so as to permit 
consideration of thermodynamic limits for M > and finite radial orbits. 

The connection with thermodynamic s is given by the equations 

Fr{ß,V,N)^~^lnZr, (40) 

where Fr is the Helmholtz free energy at radial coordinate r for the orbit of the gas 
for the metric g^ß defined by (|22]|-(|27]|. 



4 Special relativistic and Newtonian limits 

From ([T2l)-(|2T]) we see that the RaßyS vanish as r ^ «o. Thus, from (|22]l,(|23ll, (|25] |. 
(|26ll, and 

^00^^00^ /'^O, g'^~-d'\ (41) 
Po^~ \J di'^pjpk + m'^, (42) 

as expected. 

Let Zoo be the partition function as defined in (|38] | except with pQ replaced by its 
limiting value — 5-i^pjPk + »i^, as r ^ oo. The free energy for Zoo is known to be 
the Helmholtz free energy Foo{ß,V,N) for an ideal gas in special relativity. An expUcit 
calculation for Foc{ß,V,N) was given, for example, in In the next theorem, we 
show that when r = oo, so that the gas is far from the gravitational source, Fi-{ß,V,N) 
becomes the free energy of an ideal gas in special relativity. 

Theorem 2 For any finite volume V, 

]miFr{ß,V,N)=F4ß,V,N). (43) 

Proof Since f{M,r) °o as r ^ °o, the condition \x'\ < f{M,r) required by Theorem 
1 is satisfied for all sufficiently large r. The result follows easily from the Lebesgue 
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Dominated Convergence theorem for the interchange of limits and integrals because 
the integrand in (l38l l is dominated uniformly in r by an integrable function which is a 
product of terms of the form exp(— | + \p2\ + l/'sD) ■ D 

Remark 2 In Theorem 2, the Permi coordinates for the volume V relative to the orbital 
geodesic of radius r are held constant, independent of r (see Remark 1). Therefore, the 
volume of gas itself depends on r in the sense that it is held in orbit at radius r, and 
its proper volume is a function of r, converging asymptotically to its special relativistic 
volume as r ^ oo. Alternatives are also possible. The proper volume of the gas Con- 
tainer may be held fixed by choosing the limits of integration in the volume integrals 
in (l38T l to be suitable functions of r. 

Curvature in relativistic spacetimes corresponds to tidal forces in Newtonian me- 
chanics. For the purpose of identifying relativistic effects, we compute here the New- 
tonian canonical free energy and pressure of an ideal gas in a circular orbit around a 
massive body, subject only to tidal forces. 

To that end, let V be the volume of the gas orbiting a central body of mass M in a 
circular orbit of radius r. The origin of Cartesian coordinates is located at the center 
of the box with the x'-axis along the radius of the orbit. Thus, at a particular instant in 
time, the coordinates of the center of mass of the central body are (r, 0, 0). We note that 
it would be necessary to use a different coordinate System to prove Theorem 3 below, if 
we had chosen a point on the geodesic orbit, other than that corresponding to x'' = 
because of the dependence of the Newtonian limits of the Permi tetrad given by Eqs. 
Q - (0) on ^ oc x^. 

Aside from elastic collisions with the wall of the Container, the molecules, each of 
mass TO, are in free fall subject only to tidal forces. The tidal force Ft[x^ ^x^ ^x^) can be 
expressed in terms of the gravitational force Fg{x^ ,x^ ,x^). In our coordinate System, 



/v(x' ,x^,x^) 



Fg{x\x\x^)-F,{Q,Q,0), 



GMm 



(x' — r,x2,x^) 



+ 



,2 



(44) 



.((xl-r)2 + (x2)2 + (;c3)2)3/2 



The Potential energy <I> satisfying <I>(0,0,0) = and F, 



VO is given by 



<I>(x',x^x^) 



—GMm 



+ 



V(xl - r)2 + (x2)2 + (x3)2 

GMm{x^ + r) 



(45) 
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For the purpose later on of comparing $45\ with the potential energy obtained from the 
nonrelativistic limit of our metric, we approximate i45[ retaining terms to order 2 in 
the coordinates, 

GMm{2{x'f-{x^f-{x'f) 
4> . (46) 

For a volume V of an ideal gas of particles subject to an extemal field O, the 
Hamiltonian is given by, 

«=1(11 +*<"'))■ 

where x, denotes the position {x,y,z) of the /th particle. The canonical partition func- 
tion is then, 

- (2^IJ,.-P(-^^)^'^^'P' 

I ^ " i- ( / e-f'^'^^^y^^'^dxdydz] , (48) 



Klnrß J N\ \Jv 

The Helmholtz free energy, F{V,ß,N), is given by, 

^ßF{V,ß,N)=lnZiV,ß,N). (49) 

We note that it follows that the equation of State is given by a local Version of Boyle's 
Law ini, 

j3P(x)=p(x), (50) 
where p {x)dx is the probabiUty density for finding a particle at position x given by, 

NeM-ß^^ (51) 
^^''^-/^.exp(-i34>(x)Ä- ^^'^ 

Integrating both sides of (ISTT l over the volume gives {P)V =NkT, where (P) = ^ jyP{x)dx 
is the averaged pressure. 

We now consider the nonrelativistic limit of po, dSÜl l. which appears in the partition 
function We let 

A = -, (52) 

c 

and replace M by GMX^ and m by m/X in ( l3Üb . Since, in the nonrelativistic limit, we 
discard the mc^ term from the expression for the energy of each particle (recall that 
Pqc = —E), it is convenient to dehne for the /th particle the function hi{X) by 

{p^c + mc^),= (P^ + !^)=h,{X), (53) 
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and consider for the purpose of taking the nonrelativistic limit, instead of ( |38]) . the 
'reduced' partition function 

The integrability of the integrand in ( |54] | follows from the proof of Theorem 1 and 
CoroUary 1 . As before, The corresponding Helmholtz free energy, Fr{V, ß,N),is given 

by, 

-ßFi{V,ß,N)^lnZ;,{V,ß,N). (55) 
The nonrelativistic limit of hj{X) is given by 

where — pi+ P2+ p\- Eq. ( l56] l agrees with ( l46l l. 
In what follows we let G = 1 again. 

Theorem 3 The Newtonian limit of the relativistic free energy given by ( 1551 ) is the 
classical Helmholtz free energy, to second order in the space coordinates, for a gas 
subject to Newtonian tidal forces and otherwise ideal. Specifically, for any ß,N, and 
V, 

Ihn F^(y,ß,N)=F{V,ß,N), (57) 



where the right hand side of l[57\l is given by { 
Proof In view of ( l54b and ( |56] |, it suffices to prove that, 

lim / / ef'^'^^Uxdp = / / lim e^'^'^^Uxd^, (58) 

for any index /, which for convenience of notation we omit for the duration of this 
proof. Dehne s{X) = Xpq{)i). Factoring X into the Square root term in the equation 
for Po, it is easy to see that a branch of the Square root function may be selected so 
that s{X) may be extended as an analytic function of X regarded as a complex variable 
within a sufficiently small disk centered at the origin of the complex plane (depending 
on the mass term m). We note in particular that the definition of s{X) depends on the 
curvature tensor Clements given by (fT2]l-(fT7]i. With oc t = these expressions all 
become rational functions of A, except for (fT6] l which upon replacing M by MX^ is, 

3A3m3/2(^_2MA2)1/2 
^0113 — 7 T ■ .T 7s ■ (59) 

A branch of this function, analytic in a disk containing zero, exists. Thus, all the 
nonvanishing curvature Clements have extensions to analytic functions of A in a neigh- 
borhood of the origin. We denote the analytic extension of s{X) again by s{X). The 
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Taylor expansion for s centered at A = yields s(0) = —m and ^(0) = 0. It follows 
immediately that the function h{X) — (i(A) +m)/A^ is an analytic extension of h{X) 
to the domain of s{X) with a removable singularity at A = 0. We compute ^(0) = 
and, 



12mV^|^(0) = 3r^{p]+pl+plf 



+ 2Mn?r'p\ {6x\ -xj-xj) 

+ 2Mn?r'pl {2,x\ - l>x\ - Ix]) 

+ 2Mn? r^pl (8jc? - Ixl - 1x1 ) 

+ lMn?-r'^p2pT,X2XT, 

- 4piXi{p2X2 + P3X3)] 

+ 4Sm^{M)^/^r^/^p2X2Xi 
+ 48m\M)^/^,^/^pi{x]-xl) 
+ 36{Mrfm^{x\~xl) 

- 9{Mfm^{xl+xl-2x\f. (60) 
Restricting A once again to be real and positive, it now follows that 

dh 

lim— =0, (61) 

and from ( |6Ü1 ) that 

lim ^>0, (62) 

for any bounded volume V 3 {x\ ,X2,xt,), and all = p\ + p\ + p\ > K, for any suffi- 
ciently large positive constant/r(dependingony). For such dehne A = {{pi,P2,P3) G 
-.p^ < K} and denote the complement of A in by . The left band side of dSBT l 
may then be rewritten as, 

hm / / el^'^^^^dxdp + lim / / el^''^''^dxdp. (63) 

x^oJaJv i^qJa^Jv 

The integrand in the first term of (l63T l is continuous in all of its variables and is therefore 
bounded on the compact region of integration and an interval for all A ^ sufficiently 
small. From the Lebesgue dominated convergence theorem the limit may thus be inter- 
changed with the integral. In order to justify the interchange of the limit and integrals 
for the second term of ( l63T l, we observe that (A) > for all A > sufficiently small, 
and thus §jiX) > ||-(0+) = for sufficiently small A > 0. It follows that h{X) is an 
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increasing function for A > sufficiently small. The interchange of the limit and the 
Integrals for the second term in ( |63] ) now follows from the integrability of h{Xo) for 
fixed Ao > and the Lebesgue dominated convergence theorem. □ 

5 Conclusions 

There is a dearth of examples in the scientific Uterature of calculations of specific Sta- 
tistical mechanical ensembles within the broader framework of general relativity. Even 
the simplest Statistical ensemble - the canonical ensemble for an ideal gas in a finite 
volume - poses special technical problems in the context of general relativity, includ- 
ing even the existence of finite volume partition functions as finite integrals, addressed 
by the coroUary to our Theorem 1 . 

Our investigation was motivated by a simple physical question: What would be the 
gravitational effects on Statistical mechanical experiments carried out and observed in 
a spaceship in orbit around a large central mass? To test our results we found rigorous 
proofs of limiting cases: the nonrelativistic, Newtonian limit and the special relativistic 
limit for infinite radius of the circular orbit. 

Our methods apply, with minor changes, to an observer in circular orbit around a 
Reissner-Nordström charged mass. Likewise using the results of 1 1 8 1 one may consider 
the case of an observer in a general timelike equatorial orbit in the Kerr spacetime. A 
further generalization would be to include higher order terms in the spacetime variables 
in the expansion of the metric tensor in Permi coordinates. A starting point for that 
purpose are the results in ||2Ö| . 

The generalization to a gas of nonrelativistic particles with interactions, for the 
purpose of studying the effects of the background curvature, is straightforward. This 
can be achieved by following the classical formulations and including a potential en- 
ergy function of the proper distance between particles in the Boltzmann factor in the 
partition function, ( |38] |. For a more satisfactory treatment of interacting particle Sys- 
tems, it is possible that the approaches taken in the context of special relativity, such 
as in 121 and |3 |, could be applied (with second order corrections) to the circumstances 
considered in this paper, owing to the locally Minkowskian coordinates employed here. 

Returning to the ideal gas case, in Section 8 of f^, Chernikov found a covari- 
ant equilibrium single-particle Maxwell-Boltzmann distribution function expressed in 
Schwarzschild coordinates and constructed with help of an exact constant of the mo- 
tion associated with a timelike Killing vector field. That Killing field coincides with 
our {m"} on the circular geodesic orbit However, there is no exact transformation 
law in closed form from Schwarzschild to the Permi coordinates, only an approximate 
transformation law. So that approach necessarily yields approximate Statistical me- 
chanical ensembles, as in this paper. It is also unclear how the Killing field may be 
exploited to prove an appropriate Liouville theorem in Permi coordinates, which is the 
preferred coordinate System of the observer along the orbit, 
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A Invariance of the 7-form. 



Proof Consider a change of coordinates 

„0/ „1 



and 



Denote the Jacobian by, 



PO 



dx' 
dx^ ' 



j_ jf ^_ d{x'^ ,x\x^ ,X^ ,pi,P2,p-i) 

J^J{x,p) = 



(64) 



(65) 



(66) 



d{x'^,x^,x^,x'^,pi,P2,P3) ' 

Since the spacetime coordinates do not depend on momentum, it follows from a prop 
erty of block triangulär determinants that 



d{x°,x\x^,X^) d{pi,P2,P3) _ w ^ w ^ 

■^i^'P) = 3/ n I 2 3^ 37 ^=JWJ{P)- 



d{x' 



d{pi,P2,P3) 



Substituting dSÜl l for pQ in the right side of ( |65] l. we find. 



Using 



dpf, dx' p' dx^ 
dpi p^ dx~^ ' 

and the multilinearity of determinants, we may write. 



(67) 



(68) 



J{p) = -ndet 
P" 



/ 


p' 


p' 




p' 


\ 






5x1 


5x2 


5x3 






dx^ 


5? 


dx^ 


5? 








5x1 


5x2 


5x3 






5x2 


5? 


5? 


5x2 






5x0 


5x1 


5x2 


5x3 




V 


5x3 


5? 


5? 


5? 


/ 



(69) 



Next, using the relationship, p" ~ {dx" /dx^) p^ , we may rewrite the first row of the 
matrix in (|69] l so that. 



5x° 



dxß 



(70) 
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Then substituting ( ITÜl l into ( |67] ). we get, 











det 









^, (71) 



J ~ —K det 
proving the invariance of the 7-form. □ 

B Proof of Theorem 1. 

The proof of Theorem 1 depends on three Lemmas. 

Lemma 1 Let B be a 3-by-3 real Symmetrie matrix with (necessarily real) eigenvalues 
•^1 5: ^^2 5: ■^3- Thenfor any column vector p ^M? 

p^Bp>X,{p\+pl + pl) = X,\\p\\\ (72) 

where T denotes transpose and pi is the ith component of p. 

Proof Let ff be the real unitary matrix that diagonaHzes B, so that 

p^Bp = q'^ff'^Bffq 

Ai 




Az 1^ 
A3 

= Ai^f + A2^2 +■^3^3 

> Aill^lp 

= h\\ff^p\\^^Xi\\p\\\ 

where q = ff^p and where the last equation follows from the fact that unitary transfor- 
mations preserve norms. □ 

Remark 3 The quadratic form g^'^pjPk appearing in JSÜl l may be expressed in the form 
p^Bp. 

Lemma 2For any e > 0, there exists a positive function /g (M, r) defined forr> 3M > 
such that: 

(a) fs: {M, r) is a decreasing function ofM and an increasing function of r. 
Moreover, fg (M, r) °° as r ^ °° or M ^ when the other variable is heldfixed. 

(b) |1 < e, < 2e/3, \g'j\ < e/3/or / ^ j, and all eigenvalues of the 3-by- 
3 matrix whose {i,j)th component is g'^ exceed 1 — £ whenever \x^\ < f£{M,r) for 
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k=l,2,3. 



Proof Let £ > 0. From (fT2l l- (l2Tl i. it follows that by choosing | sufficiently small (for 
k= 1,2,3) depending on M and r, \Riij„,x'x"' \ < e for all /, j. Let fe (M, r) be the supre- 
mum of such bounds on The monotone and limiting behavior of /g (M,r) follows 
from the opposite monotone behavior of the Riemann curvature tensor elements and 
their zero limits (see (fT2b-(l2T]) as r ^ oo or M ^ 0. The bounds on the metric tensor el- 
ements, inverse metric tensor elements, ( |23] ), dZSl l, and ( |26] | when \x^\ < f£{M,r). The 
lower bound for the eigenvalues follows directly from the Gersgorin Circle theorem for 
the location of eigenvalues of a Square matrix flÖl . □ 

Lemma 3 For any positive real numbers a,b, and c, 

Va^+b^ + c^ > -^{a + b + c). (73) 
v3 

Proof For any x and y, x^ > 2xy. Thus, 

a^ + b^ + c"^ = i(fl^ + fo^) + i(fl^+c^) + ^(/7^ + c^) > ab + ac + bc. 
Multiplying both sides by 2/3 and rearranging terms gives, 



a^ + b^ + c^>^{a + b + cf. □ 

Proof of Theorem 1. Proof From Lemmas 1, 2, and 3 with |jc'''| < f£{M,r) for k 
1,2,3, 



\/ ig^'Pi)^ - gVpjPk ~ g^^'m^ > ^J-gVpjPk, 

> ^{^-£?{p\ + pI + pI), (74) 

> l=[i-e)(\p,\ + \p2\ + \p^\). 

Adding —g^'pi to both sides and applying Lemma 2 yields, 

> ~Y^\P1\ + \P2\ + \P3\) 

+ -L(i_e)(|^;| + 1^21 + 1/^31). (75) 

CoUecting terms, dividing by g^, and again applying Lemma 2 yields, 

PQ<-k{\Pi\ + \P2\ + \P3\), (76) 

for some positive number k provided that e < 3/(3 + 2V3). The function /(M, r) may 
be defined to be fe (M, r) for any such value of e. □ 
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